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Origin of non-keplerian motions of masers in NGC 1068 
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Abstract. Wc demonstrate that the "sub-keplerian" rotation curve of maser spots in NGC 1068 can be explained by 
the gravitational attraction of the disc orbiting the central black hole. Possible parameters matching observations 
are: black hole mass ~ (1.2 ± 0.1) x 10 7 Mq, disc outer edge > 1.3 pc, aspect ratio 3 x 10 -3 < -j| < 0.3, surface 
density E oc fl-i-OBio.i^ and digc 

mass ~ (9.4±1.6) x 10 6 M . The physical conditions required for the excitation 
of masers are fulfilled, and the outer disc would stand in a gravitationally marginally stable state. 
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1. Introduction 

Even though direct images are missing, evidences that 
Active Galactic Nuclei (AGN) host a super-massive black 
hole surrounded by a parsec-size accretion disc are com- 
pelling. In a few objects, the spatial resolution is just suf- 
ficient to unveil the outermost regions of the central en- 
gine at the sub-parsec scale thanks to the detection of 
water maser emission (Greenhill, 2002). The best exam- 
ple is NGC 4258 where masers have been associated with 
a thin disc in quasi-perfect keplerian rotation (Miyoshi 
al. 1995). In contrast, the azimuthal velocity of masers in 
NGC 1068 varies with the radius R as #-0-3i±o.02 accorc }- 
ing to Greenhill et al. (1996) who mentioned the possible 
role of the disc/torus self-gravity. To our knowledge, this 
hypothesis has never been explored, even qualitatively. We 
show in section |^ that a sub-keplerian velocity exponent 
is theoretically expected close to the outer edge of a (low 
mass) disc when considering self-gravity. In section ^, we 
report possible parameters for the disc and for the black 
hole in NGC 1068 such that the computed velocity match 
the observed motion of maser spots. A few remarks and 
possible improvements are found in the last section. 



2. Basic ideas 

2.1. Gravity, orbiting velocity and velocity exponent 

We consider the steady state system^ schematically pic- 
tured in Fig. [j], made of a black hole with mass Mbh and 
Schwarzschild radius Rs = 2 Jf bh surrounded by an axi- 

Send offprint requests to: Jean-Marc.Hure@obspm.fr 

1 In the following, quantities labeled with "bh" refer to the 
black hole, and those labeled with "disc" refer to the disc. 
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Fig. 1. Schematic view (not to scale) of the system. 



symmetrical, finite size gaseous disc with inner edge i?i n , 
outer edge R ou t, semi-thickness H, surface density £ and 
mass Af(ji sc = qM^. Assuming that gravity balances the 
centrifugal force in the disc, then the rotation velocity v$ 
of gas at cylindrical distance R from the center obeys the 
Euler equation 



-Rgi 



(>o), 



(i) 



where gn — g R h - 



-g R lsc is the radial component of the total 



(black hole + disc) gravitational acceleration. Note that 
we neglect the motion of the system around its center 
of mass (although possibly significant) and consider the 
black hole as "fixed". Further, the mass distribution in 
the disc is such that g R < at the outer edge. We can 
define the total "velocity exponent" k as 



d In Vcj) 
d\nR 



dln(-g R ) 
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(2) 



so that the azimuthal velocity varies as R k , locally. In 
particular, for a black hole only, one finds the keplerian 
value k hh 



1 



where r = \/ R 2 + z 2 is the spherical 
radius (and k hh = — | at the equatorial plane z = 0). 
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2.2. Combining a black hole and a disc: edge effects 

There is no reliable, regular expression^ for the radial field 
g^. lsc due to a disc (as a tri-dimensional object) with given 
shape and density distribution, and thus no straightfor- 
ward expression of the disc velocity exponent fc dlsc is avail- 
able. For the present system however, the relation 



<„bh 



bh 
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~disc 
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(3) 



can easily be established from Eq.(||). We thus see that 
a departure to the keplerian motion (that is k ^ k hh ) 
can arise in two opposite situations: i) when the disc self- 
gravity exceeds the gravity due to the black hole, or ii) 
when the disc self-gravity is much smaller than the cen- 
tral gravity, but |fc dlsc | ^> 1. In the extreme, the first case 
suggests that a non keplerian rotation can be produced by 
a disc only (without a black hole; i.e. q — > oo). Here, we 
shall adopt the second point of view which rather corre- 
sponds to configurations with a mass ratio q < 1. 

The point is that, depending on the mass distribution 
in a disc, large values of k dlsc can occur at the edges, 
and especially at the outer edge of interest here, even 
when q <C 1. The physical reason is that radial gravity 
is generally vanishingly small well in the middle of a disc 
because the gravitational contributions of matter located 
interior and exterior to any point essentially cancel out 
(unless very steep density gradients). Conversely, near an 
edge, there creates an imbalance between these two con- 
tributions, and the resulting field can have a quite large 
amplitude. Hence a gradient of radial self-gravity. An il- 
lustration of this "edge effect" is given in Fig. || which 
displays the radial field and velocity exponent in a sys- 
tem with mass ratio q = 0.1. In this example, we have 
adopted parameters typical for an AGN standard disc 
(e.g. Frank, King & Raine, 1992), namely black hole mass 
Mbh = 10 8 Mq, disc inner edge R m — 3i?s and outer 
edge i?out = 10 4 i?s, aspect ratio || = 0.01 x f(R) and 
surface density E cx i?~ 3 / 10 x f(R), where we have intro- 
duced the /-function in order to make the disc thickness 
and density gradually decrease to zero over a disc scale 
height typically (i.e. / ~ 1 over almost all the disc extent 
and / = outside [i?;„, Rout])- This is necessary to avoid 
a sharp discontinuity between the disc and the ambient 
medium which would not be realistic at all. Further, we 
impose that the vertical stratification of matter between 
the disc mid-plane and its surface is quadratic with the 
altitude. This conveniently mimics the Gaussian profile 
expected in a vertically isothermal disc as well as the so- 
lution of the plane Lane-Emden equation for a vertically 
self-gravitating disc (e.g. Ibahez & Sigalotti, 1984). Quite 



2 It is worth noting that any expression for gr^ ISC established 
assuming a zero-thickness disc inevitably diverges at the edges 
(Mestel, 1963) and must not be used in the disc inside and in 
the neighborhood as well. This is why we shall employ here a 
numerical tri-dimensional approach which provides everywhere 
with a regular self-gravity field. 
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Fig. 2. Example of radial field (top) and velocity ex- 
ponent k (bottom) at the equatorial plane due to a 10 s 
Mq black hole and to a low mass standard-type disc with 
outer edge located at 10 4 i?s (see text for more details). 



importantly, we point out that both the adopted prescrip- 
tion for f(R) and the assumption regarding the vertical 
stratification have no noticeable influences on the global 
results within physically acceptable limits, so that issues 
raised here are not artefacts due to discontinuity effects. 
Finally, the radial gravitational acceleration due to the 
disc is computed numerically using the accurate Poisson 
3D-solver described in Hure (2002). 

We see that, although the disc is "light" with respect 
to the black hole and central gravity fully dominates, k 
deviate markedly from the keplerian value near the outer 
edge. A general property (met here) is that the rotation is 
characterized by a "sub-keplerian velocity exponent" for 
R < R out in the sense that k > — h, whereas this is just 
the opposite for R > i? ut- 

3. The case of NGC 1068 

3.1. Origin for the non-keplerian motions 

NGC 1068 is probably one of the best studied^] among 
type-2 Seyfert galaxies. In this object with systemic veloc- 
ity of ~ 1126 km/s and distance of 15 Mpc, water maser 
emission has been detected at the parsec scale. Some emis- 
sion features would originate from the surface of a disc, 
along a diameter perpendicular to the line of sight, in the 
range 0.65 — 1.1 pc (Greenhill et al. 1996). This corre- 
sponds to a masing disc extending from ~ 9 to 15 milli- 
arcseconds (mas) in the sky. From the rotation curve of 
water maser spots, the dynamical mass enclosed within 
1.1 pc is estimated to ~ 1.5 x 10 7 Mq and the best fit for 
the velocity exponent yields fc bs. = —0.31 ± 0.02. 



See for instance the proceedings of the workshop published 
in the '97 special issue of Astrophysics and Space Science. 



Jean-Marc Hure: Origin of non-keplerian motions of masers in NGC 1068 3 



observations 
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Fig. 3. Theoretical rotation curve (bold line) expected at the surface of the disc in NGC 1068 compared to observations 
(adapted from Greenhill et al. 1996): (a) thick disc solution, (b) thin disc solution, and (c) associated velocity exponents 
(see also Tab. pi). Also shown are the contributions of the disc (plain line) and of the black hole (dashed line). 



The significant departure to Kepler law means that 
this dynamical mass may not be that of the black hole. 
Greenhill et al. (1996) suggested that disc self-gravity 
could be responsible for this non-keplerian motion. Two 
other, plausible explanations have be defended: radiation 
pressure effects (Pier & Krolik, 1992) and extra grav- 
itational attraction by a stellar cluster (Kumar, 1999). 
However, the disc self-gravity hypothesis is quite natural, 
at least because NGC 1068, as an active galaxy, hosts a 
black hole and a disc with a certain mass. It is also corrob- 



parameters 



thick disc 



thin disc 



orated by the arguments exposed in §2.2. From this point 
of view, the condition fc b s . > k hh (i-e. a sub-keplerian 
exponent) implies that the maser spots in this AGN ro- 
tate inside the disc (see Fig. |J), and not outside (which 
would mean a " super- keplerian" exponent). We can thus 
conclude that the outer edge of the masing disc is larger 
than ~ 1.1 pc. 

3.2. "Inverting" the rotation curve 

Matching observational data means that one must find 
the right magnitude of the orbiting velocity, and not only 
the velocity exponent. As v^, depends on the distribution 
of matter in space (see Eq. (Q)), we see that the prob- 
lem is to determine, if possible, which mass distribution(s) 
can cause the observed motions — a common problem in 
Astrophysics. More formally, this means the inversion of 
some relation 



Vcj, = v<f,(R, z; Mbh, Rin, -Rout, H, , 



(4) 



which obviously involves several parameters. In such a 
problem, not a single solution does exist but many, and 
finding all solutions would require the exploration of the 
whole parameter space. This task is complex and well out 



M bh (10 Y M o ) 


1.2 ±0.1 


1.2 ±0.1 


-Rout (PC) > 


1.3 (18 mas) 


> 1.5 (21 mas) 


H/R 


0.3 


« 0.3 


£ (0.65 pc) (g/cm 2 ) 


~ 470 


~ 285 


n (0.65 pc) (cm" 3 ) 


~ 10 8 


> 10 8 


s 


-1.05 ±0.1 


-1.00 ±0.05 


q (0.65 pc) 


0.49 


0.28 


q (-Rout) 


0.91 


0.65 


k 


-0.28 ±0.06 


-0.29 ±0.08 



Table 1. Parameters inferred for two "extreme" solutions 
matching observations (see also Fig||). 



of the scope of this paper. However, preliminary investiga- 
tions have revealed that i) the location of the inner edge 
has no influence provided Ri n <C -Rout, and it is set to 3Rs, 
ii) the disc thickness is not a critical parameter, provided 
the disc is geometrically thin (namely <§; 1), an d iii) 



the outer edge should stand beyond 1.1 pc (see §3.1). 

To simplify more the quest of solution, we further as- 
sume an uniform aspect ratio ^, with 0.3 as a firm upper 
limit (Gallimore et al., 1996), and we impose that the 
surface density is of the form S = T^oR s (we keep on us- 
ing the smoothing /-function; see §2.2). Besides, we set 
z = ±H, since we are interested in the velocity field at 
the disc surface where masers are assumed to be excited. 
Then, remains to be constrained mainly the central mass, 
the surface density and the location of the outer edge. 

3.3. Properties of the outer disc and central mass 

Two extreme solutions will be discussed here: a thick disc 
solution with the maximum value for aspect ratio, and a 
thin disc solution. Table [l] lists the parameters deduced 
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for the disc and for the black hole. As Fig. || shows, both 
models fit very well observations, despite the total velocity 
index k is far from being constant (but is of order of —0.3 
on average) . It turns out that i) a single value for the black 
hole mass is required, ii) S decreases roughly as J? -1 , hi) 
the disc outer edge is larger than 1.3 pc (or ~ 18 mas), 
and iv) the disc mass is close the black hole mass. Note 
that the data reported here concern the outer disc and 
must not extrapolated down to the central black hole, but 
only over one or two decade(s) in radius. 

These results deserve two comments. First, the phys- 
ical conditions required for the emission of water maser 
seem met (Neufeld, Maloney & Conger, 1994). Actually, 
at 0.65 pc, adopting £ ~ 380 g/cm 2 as a mean value (see 
Tab. |l|), the typical number density is 

s 3 x 10 7 _ 3 ._. 

(" ~ 7777 T7T\ cm > ( 5 ) 



where y, is the mean mass per particle (~ 2.3 in the molec- 
ular phase for cosmic abundances) in units of the mass 
of hydrogen tor- This is in the right density range (i.e. 

10 8 — 10 10 cm -3 ) provided the aspect ratio of the disc 
satisfies 3 x 10~ 3 < |j < 0.3. Further, for such densities 
and for temperatures appropriate for molecules to survive 
(i.e. < 3000 K), support by radiation pressure must be 
significant, even dominant (whatever the solution). This 
conclusion is akin to the analysis by Pier & Krolik (1992). 
Assuming a vertically self-gravitating disc in hydrostatic 
equilibrium, the mid-plane temperature is 



V 4^ 



1/4 



2000 K, 



(6) 



where a is the Stephan constant and c is the speed of 
light. For a typical grain opacity of 1 cm 2 /g, the vertical 
optical depth is r > 200, implying a temperature at the 
disc surface of ~ T^d/r 1 / 4 ~ 600 K, as required (Neufeld, 
Maloney & Conger, 1994). 

The second remark concerns the gravitational stability 
of the disc. For both solutions, the Toomre Q-parameter 
is of the order of unity. Actually, again at 0.65 pc, we have 



D KCs ~ V^nid . l 8crH 9 r-s 
V ~ ttG£ ~ nGR V3c£3 ~^ V) 

where k ss f2 is the epicyclic frequency and c s is the sound 
speed. It means that the outer disc is in a marginally sta- 
ble state, and does not resemble a standard disc (Kumar, 
1999). This property is predicted by the standard model 
at the parsec scale, whatever the accretion rate in the in- 
ner disc. Note also that the fact that our solutions have 
£ oc R^ 1 and v^ lsc ~ est (see Fig. ||) is fully consistent 
with what we expect in an (radially) isothermal Q — est 
disc (e.g. Mestel, 1963). 



4. Concluding remarks 

We have reminded that a non-keplerian motion can arise 
from the combination of the gravitational attraction of a 



black hole and a disc. We propose that the "sub-keplerian" 
profile traced by water masers in NGC 1068 results from 
such an effect, a possibility outlined already in Grecnhill 
et al. (1996). We have reported two simple physical solu- 
tions producing approximately the observed mean velocity 
exponent, with, accordingly, the parameters for the black 
hole and for the outer disc in this active galaxy. As out- 
lined, these solutions are not unique; for instance, it is 
possible to fit the feature near 12 mas (see Fig. ||), but £ 
is no more a power law. From this point of view, it would 
be worthwhile to explore the whole parameter space. 

The analysis discussed here can be applied to other 
AGN. Note however that the system considered here 
is very simple and would need refinements. In partic- 
ular, there is no physical model for the disc and the 
present approach works with the most minimal condi- 
tions. For instance, Eq.(|l|) should be supplemented with 
non-gravitational terms (pressure gradient, advection, vis- 
cosity, etc.). Further, the gravitational attraction, possi- 
bly significant, of other components of the AGN (molecu- 
lar torus, stellar cluster, galactic bulge, etc.) is not taken 
into account. Significant deviations to axi-symmetry (like 
a warp) could also change our results, quantitatively, and 
would make the problem much more tricky. With these 
restrictions, we have shown that it is feasible to "invert" 
the rotation curve in order to constrain the mass distribu- 
tion in the central engine at the parsec scale. With more 
sophisticated or self-consistent investigations (e.g. Lodato 
& Bertin, 2002), one should be able to connect these outer 
solutions to the innermost ones where accretion theories 
are very uncertain. 
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